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We review the theoretical basis and understanding of electronic interactions in graphene Landau 
levels, in the limit of strong correlations. This limit occurs when inter-Landau-level excitations may 
be omitted because they belong to a high-energy sector, whereas the low-energy excitations only 
involve the same level, such that the kinetic energy (of the Landau level) is an unimportant constant. 
Two prominent effects emerge in this limit of strong electronic correlations: generalised quantum Hall 
ferromagnetic states that profit from the approximate four-fold spin- valley degeneracy of graphene's 
Landau levels and the fractional quantum Hall effect. Here, we discuss these effects in the framework 
of an SU(4)-symmetric theory, in comparison with available experimental observations. 
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I. INTRODUCTION 



The theory of non-interacting masslcss Dirac fermions 
in two spatial dimensions provides the framework which 
allows for the understanding of most of graphene's low- 
energy electronic propertiesii At first sight, this may 
seem astonishing because the Coulomb interaction is 
strictly speaking of an intermediate strength; indeed, 
within a typical Coulomb-gas argument, one compares 
the average interaction energy e^kp/e, at the character- 
istic length ~ kp^, in terms of the Fermi wave vector 
kp and the dielectric constant e of the environment sur- 
rounding the graphene sheet, to the kinetic energy Tivpkp 
at the same length scale. The ratio between these ener- 



gies yields the coupling constant ac = e /hevp ~ 2/e, 
which is reminiscent of the fine-structure constant in 
quantum electrodynamics if one replaces the speed of 
light c by the Fermi velocity vp, i.e. the characteristic ve- 
locity of the electrons in a material. Because w^? ~ c/300 
in graphene, the graphene fine- structure constant ac is 
roughly 300 times larger than that (a ~ 1/137) of quan- 
tum electrodynamics. 

In view of the rather large coupling constant, one might 
expect to see correlation effects in graphene, which hap- 
pen though to be sparsci^ Indeed, electronic instabili- 
ties are not only triggered by the (bare) coupling con- 
stant, but one also needs to take into account the den- 
sity of states at the Fermi level in the discussion of such 
instabilities j^i^ As a consequence of the linearity and the 
two-dimensional (2D) character of graphene electrons, 
the density of states, however, vanishes linearly with the 
Fermi energy when approaching the limit of undoped (in- 
trinsic) graphene, such that electronic instabilities are 
suppressed. In the search of prominent correlation ef- 
fects, one should therefore investigate situations in which 
the density of states in graphene is enhanced. A (loga- 
rithmically) diverging density of states is typically en- 
countered at van Hove singularities due to saddle points 
in the band dispersion. Van Hove singularities occur at 
extremely high energies (^ 3 eV) in monolayer graphene 



and arc thus inaccessible with the help of field-effect 
doping. In contrast to monolayer graphene, they occur 
at rather small energies in AB-stackcd bilayer graphene 
(^ 3 mcV),— which arc only resolved at electronic den- 
sities below 10^^ cm~^. A promising system in this re- 
spect is twisted bilayer graphene, where van Hove singu- 
larities at intermediate energies (^ 10. ..100 meV) have 
been observed.— Such twists naturally occur in epitaxial 
graphene grown on the carbon face of the SiC crystalii 

An alternative means of inducing a large density of 
states in graphene, and thus of increasing the role of elec- 
tronic correlations, is to expose the sample to a strong 
perpendicular magnetic field B. In this case the elec- 
tronic energy is quantised into highly degenerate Lan- 
dau levels (LLs) at discrete energies En, around which 
the density of states p{E) is strongly peaked, p{E) = 
guB J2n /(^ ~ En), where ub = eB/h is the flux den- 
sity measured in units of the flux quantum h/e and g 
takes into account internal degrees of freedom, such as 
the four-fold spin-valley degeneracy {g = 4) in graphene. 
The functions f{E — En), which are normalised to one, 
/ dE f{E~ En) — 1, become delta functions in the clean 
limit, which we assume in the theoretical discussion here. 
Each LL may thus be viewed, in this limit, as an infinitely 
flat energy band the density of states of which grows lin- 
early with the magnetic fleld. 

In this article, we review some effects due to 
the magnetic-fleld induced electronic correlations in 
graphene, in comparison with the perhaps better-known 
2D electron gas in semiconductor heterostructures. The 
probably most prominent one is the fractional quantum 
Hall effect (FQHE), which has recently been observed ex- 
perimentally in the two-terminal^ and the four-terminal 
configuration] ^°d^ In contrast to the FQHE in GaAs 
heterostructures, the graphene FQHE reflects a four- 
component structurei^ii^ that is inherited from the four- 
fold spin- valley degeneracy and that goes along with par- 
ticular magnetic properties described in the framework of 
SU(4) quantum Hall ferromagnetism.— The latter is also 
relevant in the discussion of interaction-induced integer 
quantum Hall effects (IQHE) at integer fllling factors v 
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that do not belong to the "magic" series v = n^ijuB = 
±2, ±6, ±10, in terms of the carrier density rig/. 

The article is organised as follows. In Sec. HIl we re- 
view some of the experimental findings from the observa- 
tions of the IQHE in 2005 to the very recent ones of the 
FQHE in the four-terminal geometry, in 2010. After an 
introduction to the theoretical basics of graphene LLs in 
Sec. Uni wc discuss the SU(4) -spin- valley quantum Hall 
fcrromagnctism in Sec. IIII Al and the SU(4) theory of the 
FQHE in Sec. HHB] 

II. EXPERIMENTAL SITUATION 

A. Relativistic integer quantum Hall effect 

A milestone experiment in graphene research was the 
observation in 2005 of a particular ~ relativistic - IQHE 
in graphene, when changing either the electronic den- 
sity via the electric-field effect at a fixed magnetic field 
or when varying the field at fixed electronic density i^^'^^ 
The samples used in these magnetotransport measure- 
ments were obtained with the help of the exfoliation 
technique^ and the effect has later (in 2009) been con- 
firmed in epitaxial graphene samples which have also 
been proven to be promising for metrological means be- 
cause of a high-precision (with an error bar on the order 
of 10^) Hall-resistance quantisation^^ Whereas the effect 
has the same signature - a plateau in the Hall resistance 
accompanied by a vanishing longitudinal resistance - as 
that in conventional 2D electron systems, it occurs at 
unusual filling factors, 

t/ = ±2(2n+l) = ±2,±6,±10, (1) 

and refiects the relativistic nature of the charge carriers 
in graphene. Indeed, the two possible signs ± reflect the 
presence of a conduction band (-1- for "particles") that 
touches the valence band (— for "anti-particles" on the 
hole-doped side) . The filling- factor steps in units of four 
between successive plateaus may easily be understood 
as a consequence of the four-fold spin- valley degeneracy, 
which was not resolved in these first experiments and 
that yields four copies of each LL. The offset of ±2 in 
the plateau series ((T]) is a consequence of relativistic LL 
quantisation that yields a LL spectrum 

Ex,,, = h^^, (2) 

where Ib = ^jTi/eB ~ 26/^B[T] nm is the magnetic 
length, and the integers n label the levels in the con- 
duction band (A = -I-) or in the valence band (A = — ). 
The most prominent feature of the level spectrum 
apart from its square-root dispersion with the magnetic 
field and with n, is the presence of a zero-energy LL for 
n = 0. It is this level that is responsible for the offset ±2 
in the series ([T]) because it is only half-filled at zero dop- 
ing, J/ = 0. This means that the condition for the IQHE, 



namely a set of completely filled LLs with a topmost filled 
level separated by a gap from the lowest unoccupied LL, 
is not fulfilled at = 0, but only at = 2 (for electron 
doping) or = — 2 (for hole doping), as a consequence of 
the four-fold spin- valley degeneracy of the n = LL. 



B. Additional plateaus at integer fillings 

In 2006, one year after the discovery of the graphene 
IQHE, novel high-field plateaus have been observed at 
ly = 0, ±1 and ±4 that do not belong to the series ([T|)i^ 
These additional states indicate that the spin- valley de- 
generacy in the n = LL is completely lifted, whereas 
in n = 1 it is only partially lifted - if it were fully lifted 
in the latter case, one would also expect an IQHE at 
1/ = ±3 and ±5. The explanations which have been 
given for the spin-valley degeneracy lifting fall into two 
classes: (1) extrinsic or (2) intrinsic, i.e. interaction- 
induced, effects. The simplest extrinsic effect is certainly 
the Zeeman effect that would lift the spin degeneracy, 
such that each four-fold degenerate LL is split into two 
(valley-degenerate) spin branches separated by an en- 
ergy scale of Az — 1.2B[T] K. A more subtle extrin- 
sic effect, as a consequence of electron-phonon coupling, 
is capable of lifting the valley degeneracy in the zero- 
energy LL in form of the generation of a mass gap in 
the level spectrumi^"— The coupling to an out-of-plane 
phonon can yield a Peierls-type distortion and can thus 
break the inversion symmetry of the latticci^i More re- 
cently an inplane Kckulc distortion has been investigated 
that couples the two different valleys i^i^ In contrast to 
the out-of-plane distortion, the latter mechanism yields 
a mass term (and thus a valley splitting) that does not 
depend on the coupling to the substrate and that has 
been evaluated to be roughly A^ek — 2B[T] K. Notice 
that the linear B-field dependence simply reflects the fact 
that the coupling is proportional to the density of states 
p{E) cx B, which scales linearly with the magnetic field, 
as mentioned in the introduction. 

The second class of degeneracy-lifting effects contains 
mechanisms that are triggered by the Coulomb inter- 
action between the electrons and that are discussed in 
more detail in Sec. IIII Al One effect is the so-called 
magnetic catalysis, which has been investigated before 
the discovery of the IQHE in graphene, in the con- 
text of Dirac fermionsi^ The mechanism consists of a 
mass-gap generation that yields the same level spectrum 
(and thus a valley-degeneracy lifting) as that due to the 
above-mentioned Peierls-type distortions, but it is dy- 
namically generated by the electron-electron interactions 
themselves 1^1^ The mass term plays the role of an order 
parameter that has been identified with exciton conden- 
sation. Independently, quantum-Hall ferromagnetism, 
both in the spin and in the valley channel, has been pro- 
posed in 2006 as a possible route to understanding the 
additional platcausj^"— Quantum-Hall ferromagnetism 
is an exchange effect, where the electron-electron interac- 
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tion is minimised by the formation of a maximally anti- 
symmetric orbital wave function, accompanied by a max- 
imally symmetric valley-spin part. The effect is partic- 
ularly efficient in LLs because the latter may be viewed 
as infinitely fiat bands - the polarisation of the spin and 
the valley pseudospin is therefore not accompanied by 
a cost in kinetic energy. Whereas quantum-Hall ferro- 
magnetism is capable of generating a transport gap, and 
thus an IQHE, at all integer filling factors that do not 
belong to the series ([T]))^ a mass gap resulting from a 
lattice distortion or magnetic catalysis can only lift the 
valley degeneracy in the zero-energy LL n = 0. However 
in both cases, magnetic catalysis and quantum-Hall fer- 
romagnetism, the gap scales with the typical interaction 
energy e^/els oc a/B- 



energy 


value for arbitrary B 


for S = 25 T 


Az 


1.2B[T] K 


30 K 


Akek 


2B[T] K 


50 K 


e^/els (vacuum) 


139 ^B[T] K 


694 K 


e^/elB (on SiOa) 


104yB[T] K 


521 K 


e'^/elB (on h-BN) 


109 ^B[T] K 


543 K 


e'^/els (on SiC) 


ny^B[T] K 


355 K 



TABLE I: Tab. I. Energy scales for spin-valley degeneracy lifting 
in graphene LLs. The first two lines show the energy scales associ- 
ated with extrinsic effects (Zeeman effect and Kekule-type lattice 
distortion, Az and Afcgj., respectively), which are proportional to 
B. Below are shown the interaction-energy scales (oc \/B), differ- 
ent substrates taking into account both the dielectric constant of 
the substrate and RPA contributions from inter-band processes. 



The discussed energy scales are summarised in the ta- 
ble above. The interaction energy scales depend on the 
dielectric constant e of the environment, which consists in 
a typical experimental situation on the substrate [Si02, 
hexa boron nitride (h-BN) or SiC for epitaxial graphene] 
on one side and air (vacuum) on the other one. The third 
line (vacuum) indicates the energy scale for freestanding 
graphene. A recent theoretical study proposes to engi- 
neer the short-range part of interaction potential via a 
partial screening with a dielectric medium at a finite dis- 
tance from the graphene sheet Notice that in all cases, 
we have taken into account screening due to the com- 
pletely filled valence band, which yields e = e(l -t-Trac/S) 
within the random-phase approximatioui^ One notices 
from these energy scales that in all cases the Coulomb in- 
teraction sets the leading energy scale and should thus be 
considered first in the discussion of the spin- valley degen- 
eracy lifting, whereas extrinsic effects are subordinate. 
As it is discussed below in Sec. IHI Al the extrinsic ef- 
fects are cooperative with quantum-Hall ferromagnetism 
in the sense that they orient the interaction-induced spin- 
valley magnetisation in a particular direction. 



C. Fractional quantum Hall effect in graphene 

In the previous subsection, we have argued that the ap- 
pearance of IQHE plateaus, which do not match the se- 
ries ([T]), could in principle be understood without invok- 
ing electron-electron interactions - extrinsic effects could 
be responsible for these plateaus, although this is unlikely 
in view of the different energy scales involved. Clear ev- 
idence for interaction-induced phases in graphene LLs 
has been found in 2009 with the first observations of the 
FQHE at 1/ = 1/3 in suspended graphenci^ One no- 
tices that it took roughly twice as long in graphene be- 
tween the observation of the IQHE and the FQHE as 
compared to conventional 2D electron systems, where 
the IQHE was observed in 1980^ whereas the FQHE 
was discovered in 1982,^ The necessary mobility increase 
of graphene samples, which is required for the observa- 
tion of the FQHE, could already be achieved in 2008 in 
current-annealed suspended samples, where mobilities in 
the 100 000 cm^/Vs range have been reportedi^^ How- 
ever, it turned out to be an experimental challenge to 
obtain samples with working and sufficiently separated 
electronic contacts. 

The above-mentioned transport measurements, which 
revealed the FQHE, were indeed performed in the two- 
terminal configuration, where the same contacts used as 
source and drain serve for the resistance measurement. It 
is therefore not possible to perform simultaneously a mea- 
surement of the longitudinal and the Hall resistance, but 
both are superposed, and sophisticated conformal map- 
pings are necessary to separate the two componentsj^ 
However, because of the vanishing longitudinal compo- 
nent in the case of the FQHE, the two-terminal resistance 
is then determined by the quantised "Hall" resistance and 
therefore reveals the characteristic plateau. 

These first observations have since been confirmed in 
the more robust four-terminal configuration, which al- 
lows for a simultaneous measurement of and thus a clear 
distinction between the longitudinal and the Hall resis- 
tances. Two experiments were reported in 2010, one 
in a suspended graphene samplers and another one in 
graphene on an h-BN substrateii that allows for a mo- 
bility increase upon current annealing that is in the 
same range as (though somewhat lower than) that in 
suspended graphenci^ In Ref. [13 the activation gap 
could be determined and agrees rather well with that 
0.05...0.1eVeZB ~ l...\Ay^B[r\ K one expects^i^s for 
the polarised Laughlin state at v = 1/3.— In graphene 
on a h-BN substrate, most members of the 1/3 FQHE 
family va n = Q {&i v — ±1/3, ±2/3, and ±4/3) and all 
in n = 1 (at v = ±7/3, ±8/3, ±10/3, and ±11/3) could 
be resolved to great accuracy— Whereas the absence or 
relative weakness of the ±5/3 member of the 1/3- family 
in the zero-energy LL remains to be understood, it clearly 
corroborates the approximate SU(4) symmetry due to the 
four-fold spin-valley degeneracy underlying the FQHE in 
graphene fi^iii^ as discussed in more detail from the the- 
oretical point of view in Sec. IHI Bl Another relevant 
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finding in graphene on an h-BN substrate is that of addi- 
tional IQHE plateaus at z/ = ±3 and ±5^^ which indicate 
a full spin- valley degeneracy lifting not only in n = but 
also in n = 1. Whereas these additional plateaus can- 
not be explained in the framework of a mass-gap genera- 
tion, either by a lattice distortion or interaction-induced 
magnetic catalysis, it is expected from the formation of 
quantum-Hall ferromagnetic states, as mentioned above. 

III. THEORETICAL UNDERSTANDING 

From the theoretical point of view, both the FQHE 
and quantum-Hall fcrromagnetism require three essen- 
tial ingredients: (1) infinitely flat and highly-degenerate 
energy bands (here in the form of LLs the degeneracy 
of which is characterised by the flux density ns); (2) the 
Aharonov-Bohm effect, which yields a geometric phase to 
paths in the 2D plane and that induces the so-called mag- 
netic translation group; and (3) sufficiently short-range 
interactions)^ 

In order to make transparent the above statements, we 
consider a single LL that is sufficiently well separated in 
energy from its adjacent levels. This condition is fulfilled 
when the LL spacing A„ = \/2h{vF /lB)Wn + 1 — ~ 
hvp/V^nlB ^ 200\/S[T]/n K is larger than the im- 
purity broadening of the levels. Furthermore, because 
this energy scale is much larger than the extrinsic spin- 
valley symmetry-breaking effects (see Tab. I), we may 
consider each LL as four-fold degenerate, in addition to 
the orbital degeneracy given by the flux density n b ■ The 
electron-electron interactions may then be separated into 
a low-energy and a high-energy part. The latter consists 
of interaction-induced inter-LL transitions at the char- 
acteristic energy scale A„, whereas the low-energy part 
consists of intra-LL excitations, in which case the kinetic 
energy (set by the scale Hvf/Ib) effectively drops out of 
the problem. The low-energy interaction Hamiltonian (in 
reciprocal space) may then be written as 

Hn = ^Yl ^il)Pnh(l)Pn{q), (3) 

q 

where v{q) — lixe^ jeq is the Fourier-transformed 
Coulomb interaction potential, and the Fourier compo- 
nents of the density operator /9n(q) take into account 
only states within the ji-th LL. Apart from a form fac- 
tor Tn{q) that takes into account the overlap between 
electronic wave functions in the n-th LL and that may 
be absorbed into an effective interaction potential)^ the 
(projected) density operator 

is a sum of the one-particle density operators Pj(q) = 
exp(— iq • Rj) for each of the N electrons in the n-th LL. 
Here, Hj = {Xj,Yj) is the operator that describes the 



centre of the cyclotron motion (called guiding centre) of 
the j-th electron. Because it is a constant of motion, 
it docs not connect states in different LLs, in agreement 
with the construction of the model ([3]). Furthermore, the 
components of Rj do not commute, 

[Xj,Yj'] = ilsSj,]', (5) 

which is a manifestation of the above-mentioned 
Aharonov-Bohm effect. Indeed, one sees from the com- 
mutation relations that Xj and Yj are conjugate vari- 
ables, such that Xj generates a translation in the —y- 
direction, whereas Yj generates one in the a;-direction. 
As a consequence the electron, when moving on a closed 
path around an area S, picks up an Aharonov-Bohm 
phase ip = ~ 27r0/(/)o, where 4> = BYi is the flux 

in the area S. The commutation relations ([S]) further- 
more induce the commutation relations 

[p-(q), p-(k)] = 2^ sin ^.ih^^il)^ + k) (g) 

for the projected density operators, such that their 
Hcisenbcrg equations of motion ih,p{q) = [p{q),H] be- 
come highly non-linear. 

Another consequence of the commutation relations , 
which allow for the introduction of harmonic-oscillator 
ladder operators bj ~ {Xj + iYj ) / y/2l b and 6^ = {Xj — 

iYj)/\/2lB, with [bj,b'j,] = Sjj', is the representation of 
n = 0-LL states in terms of analytic functions 

0„(z„z;) ^ ^ zf e-l^^l'/^ (7) 

where we have defined the complex position Zj = {xj — 
iyj)/lB of the j-th particle in the 2D plane and where 
we have omitted the normalisation constant in the ex- 
pressions. Whereas this statement is, strictly speaking, 
true only in the n = LL, one may nevertheless use 
a representation of states in other LLs in terms of ana- 
lytic functions if one interprets Vn{q) = ''-^{l)\-Fn{<i)\'^ as 
an effective interaction potential that mimics the 7i-th 
LL while considering the projected density p(q) as one 
of 71 ~ 0. This assumption is justified by the fact that 
the commutation relations (jB]) do not explicitly depend 
on the LL index. 

To summarise this theoretical introduction of the 
model, one notices the following important points. 

• The above arguments arc valid for any type of Lan- 
dau quantisation and not restricted to the rcla- 
tivistic one in graphene. From this point of view, 
graphene and its FQHE are not so different from 
the FQHE in semiconductor heterostrucures. 

• The specificity of graphene and relativistic LL 
quantisation is revealed rather in the form factors 
Tn{q), which take into account the wave function 
overlaps in a particular LL n and that happen to be 
different from that in non-rclativistic 2D electron 
systems, as a consequence of the spinorial struc- 
ture of the electronic wave functions in graphene. 
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wave vector 



FIG. 1: Spin-valley waves (Goldstone modes) in graphene LLs. 
(a) For one filled spin- valley subbranch (k = 1), there are three 
Goldstone modes connecting the filled subbranch with the empty 
ones, (h) For two filled subbranches (k = 2), four Goldstone modes 
connect the completely filled subbranches with the empty ones, (c) 
Dispersion relation of the spin- valley waves (continuous line). At 
qIb <C 1, the energy disperses as Eq oc q^, as one expects for 
magnons. At qlg ^ 1, the dispersion relation saturates and may 
be approximated by the energy to create a well-separated electron- 
hole pair, with a Coulomb attraction between the electron and the 
hole, Eq ~ [2^/^- l/qlB]e^/dB (dashed line). 



• Another specificity of graphene is its four- fold spin- 
valley degeneracy. The Coulomb interaction natu- 
rally commutes with the electronic spin, whereas 
this is a priori not the case for the valley pseu- 
dospin - indeed, the two different valleys can be 
coupled by scattering due to short-range compo- 
nents of the Coulomb potential. However, one 
may show that this intcr-valley coupling is sup- 
pressed by a factor oi q/Ib ^ 0.005y^i?[T] because 
of the reciprocal-space distance ^ between the 
two valleys;22i2^ such that the Coulomb interaction 
in graphene LLs may be viewed as approximately 
SU(4)-symmetric. 

• A general A^-particle wave function in graphene LLs 
must therefore be described by an analytic polyno- 
mial in all particle coordinates z^if\ where the su- 
perscript (j) indicates one of the four-spin valley 
components (/v,t), (K^i), (if',t), or 



A. Interaction-induced integer quantum Hall effect 

A first manifestation, the understanding of which turns 
out to be instructive also for the FQHE, of Coulomb 
interactions in graphene LLs is the formation of SU(4) 
quantum-Hall ferromagnetic states at integer filling fac- 
tors that do not correspond to completely filled LLs de- 



scribed by the series ([T]). As already mentioned above, 
a repulsive interaction such as the Coulomb potential 
favours orbital wave functions with nodes when two par- 
ticles approach each other. These nodes are naturally 
built in in the usual Slater determinants for k completely 
filled spin- valley branches of the last (partially) occupied 
LL, 

K Nb-1 

|FM) = n n 4,, I vac), (8) 

j — l m— 

where c]^j ^ creates an electron in the component j in the 
state corresponding to the LL wave function z™ and |vac) 
denotes the fcrmion vacuum. For illustration, we con- 
sider the state ([5]) in the zero-energy LL n = and omit 
the LL index at the fermion operators, keeping in mind 
that the generalisation to other LLs is straight-forward, 
as discussed above. 

Naturally, the state (|8]) is the ground state of a model 
with no interactions if the filled states have a lower one- 
particle energy than the empty ones, for example in the 
presence of a Zeeman effect. Here, however, we consider 
these effects to be absent, and we do therefore not specify 
whether the occupied branches are particular spin or val- 
ley states. The state ([5]) therefore breaks the SU(4) spin- 
valley symmetry, which is respected by the interaction 
model. As a consequence of this symmetry breaking, the 
quantum-Hall ferromagnet ([5|) has low-energy excitations 
in form of Goldstone modes the energy of which vanishes 
in the small-wave-vector limit. These Goldstone modes 
are spin-valley waves that connect the different possible 
ground states, which are obtained by relabeling the occu- 
pied subbranches. They are depicted in Fig. [T]for k = 1 
(a) and k = 2 (b) completely filled subbranches ~ the case 
of K = 3 filled subbranches is particle-hole-symmetric to 
K = 1. For K = 1, the Goldstone modes are three-fold 
degenerate [see Fig. [TJa)] whereas for k = 2 there are 
four different types [Fig. [Hb)]. 



1. Spin-valley waves 

The different spin-valley waves may be obtained by 
application of the operator 

p.,(q)= (9) 

on the state |FM), where R is the one-particle oper- 
ator associated with the guiding centre, as discussed 
above, and j denotes an occupied spin-valley compo- 
nent, whereas i corresponds to an unoccupied one. The 
state pij (q) |FM) may also be viewed as a superposition of 
particle-hole excitations, where q is the wave vector of the 
excitation. As a consequence of the magnetic translation 
algebra, generated by the commutation relations ([S]), this 
wave vector is proportional to the distance AR = R — R' 
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between the guiding centre R of the electron and that R' 
of the hole, 

q = AR X e,//|. (10) 

The energy spectrum of the spin-valley waves may be 
obtained by evaluating the Hamiltonian ([3]) in the state 
Pij(q)|FM), and one obtains 



= (FM|p,j(-q)i/„Py(q)-i/„|FM) 

^' 2 ( Qxky ^ kxQy 



2^z;„(k) sin^ 



or explicitly, in the zero-energy LL n ~ q ^28, 29,41, 42 



1 _ e-«''-/4/o 



q ii 



(11) 



(12) 



which is plotted in Fig. [ijc). In the last expression, 
which is independent of the number k of filled spin- valley 
branches and which thus indicates that the spin-valley 
waves are degenerate for all values of q, Jq {x) is a modi- 
fied Bessel function. The limits of the dispersion (fT2|) are 
transparent; for small values of the wave vector qls <C 1, 
one obtains the usual dispersion expected for spin- 
wave-type modes. 



P — ^"^2/2 



in terms of the spin stiffness 



-^Vz;o(k)|k|^Z| = ^=^. 
uv yi I a 16^2^ e^s 

k 



(13) 



(14) 



In the opposite limit, the dispersion may be under- 
stood in terms of the energy of a spatially well-separated 
electron-hole pair. The energy to add an electron (or a 
hole) to the state ([5]) is just given by the exchange energy. 
Ex = ■\/7r/8(e^/e/B). The value at which the dispersion 
([T2|) saturates is indeed twice the exchange energy. Fur- 
thermore, the electron and the hole with opposite charge 
interact via the Coulomb attraction 



(15) 



elARI 



eqll 



as a consequence of the connection ([T0| between the wave 
vector and the distance between the guiding centre of the 
electron and that of the hole. As depicted by the dashed 
line in Fig. [Ijc), the spin- valley- wave dispersion is well 
approximated by the sum of these two terms. 



between the electrons without the need of explicit spin- 
valley symmetry-breaking terms, such as the Zeeman ef- 
fect. The state is stable because the dispersion of the 
collective excitations is gapped for any non-zero value of 
the wave vector, and the addition of an electron (or a 
hole) is associated with an energy cost given by the ex- 
change energy, which is much larger than the external 
symmetry-breaking fields (see Table 1). The role of such 
external terms is then reduced to a simple orientation of 
the interaction-induced spin-valley magnetisation, simi- 
larly to a usual (spin) ferromagnet placed into a magnetic 
field that orients its magnetisation in the direction of the 
field. 

In the following paragraph, we argue that there 
are lower-energy elementary excitations, in form of 
skyrmions, than such additional electron with a flipped 
spin or valley-pseudospin. However, their energy is also 
determined by the interaction-energy scale e^/eZs, such 
that the overall picture remains unaltered. 



2. SU(4) skyrmions 

In the previous paragraph, we have considered the 
elementary excitation to be a simple additional elec- 
tron (or hole) that is added into an unoccupied spin- 
valley component in the quantum-Hall ferromagnet ([8|). 
Its energy is then simply given by the exchange energy 



However, it turns out to be energet- 



Ex = 

ically favourable for this additional particle to be dressed 
by a local deformation of the SU(4)-ferromagnetic back- 
ground, such as to lower the energy cost due to the op- 
posite spin orientation of the particle with respect to 
the background. This dressed particle is called skyrmion 
and carries a topological charge in addition to its electric 
one. For a simple SU(2) spin ferromagnet, this topolog- 
ical charge may be viewed as the number of times the 
(normalised) local magnetisation wraps, when exploring 
the 2D plane, the Bloch sphere the points of which rep- 
resent the orientation of the magnetisation. The SU(4) 
spin- valley case is more complicated and requires the in- 
troduction of two additional Bloch spheres (one for the 
valley-pseudospin and one for the two angles that de- 
scribe the entanglement between the spin and the valley- 
pseudospin))^ but the picture is essentially the same. 

The topological charge Qtop, which is a positive or 
negative integer, determines the energy of the skyrmion 
excitationJ^i^ 



E,k = 47rps|(3 



topi 



(17) 



E, 



qlB»l 



1 



8 qh 



e 

els 



(16) 



To summarise the picture of SU(4) quantum Hall fer- 
romagnetism and the associated spin- valley-wave modes, 
we first mention that the polarised state ([5]) may be ob- 
tained simply as a consequence of the Coulomb repulsion 



in terms of the spin stiffness ([T4| . One thus notices that 
the energy to create a skyrmion with charge Qtop = ±1 
is half of that to create a simple (undressed) electron 
in an unoccupied spin-valley component. Dressing this 
additional electron by a topological spin-valley texture 
therefore lowers the energy in n = by a factor of 2. 
The energy gain is less in the LLs n ^ 0, but it remains 
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positive for n = 1 and n = 2, whereas in even higher 
LLs it is no longer energetically favourable to dress the 
additional charge by creating skyrmionsj22i22 

We finally mention that the skyrmion is generically 
larger in size than an undressed electronic excitation 
Ib)- The size of the skyrmion is indeed determined 
by a competition between the Coulomb (exchange) in- 
teraction, which favours large skyrmions to maintain lo- 
cally the ferromagnetic order, and external symmetry- 
breaking terms that, even if they are small, have a ten- 
dency to lower the number of reversed spins or valley- 
pseudospins and thus to lower the skyrmion size. Indeed, 
the skyrmion radius scales as^i^ 

^-]l^lB, (18) 

where A represents a generic spin-valley symmetry- 
breaking terms, such as the Zceman effect or that arising 
from a spontaneous lattice distortion discussed in Sec. 

EH 



B. SU(4) fractional quantum Hall effect 

In the previous section, we have argued that electron- 
electron interactions are responsible for the formation of 
additional plateaus at integer filling factors that do not 
correspond to the series (P), as a consequence of the for- 
mation of maximally polarised quantum-Hall ferromag- 
netic states. These considerations turn out to be helpful 
also in the understanding of the four-component FQHE. 
If we consider, e.g., the SU(4) ferromagnetic state at 
I' = ± 1 , its orbital wave function may be written in terms 
of the completely anti-symmetric Slater determinant 

N ^ 

cl>{{zk}) ^Y[{zk- zi)e~^>^-^\''\"/\ (19) 

k<l 

in terms of the complex coordinates Zk of the electron in 
units of Ib, regardless of the spin- valley component they 
belong to. As a consequence of the anti-symmetry of this 
orbital wave function, the associated spin-valley wave 
function must be completely symmetric, i.e. precisely 
ferromagnetic, such as to fulfil the anti-symmetry re- 
quirement for fermionic A'^-particle wave functions. Wave 
function ([T9l) is the simplest example of Laughlin's wave 
functior^ 

N ^ 

<l^t {W}) ^X{{z,~ zir e (20) 

k<l 

which describes FQHE states at filling factors v = 1/m. 
Indeed, a power counting of the terms in the polynomial 
indicates that the largest power of an arbitrarily chosen 
particle component Zk is M = m{N — 1). As we have 
already mentioned in the first part of this section, this 



power is delimited by the number of flux quanta thread- 
ing the 2D system, such that M = Nb — 1, and one 
obtains, in the thermodynamic limit, the relation 

Nb ub 1 , . 

= -irr = — = (21 

i.e. the exponent m in Laughlin's wave function deter- 
mines the filling factor. For odd values of m - remember 
from the previous discussion that m must be an integer 
to match the analyticity condition for wave functions in 
the LL rt = - the same symmetry arguments apply as 
for the wave function ()19p . It is a fully anti-symmetric or- 
bital wave function, and the spin-valley part must there- 
fore be completely symmetric, such that Laughlin's wave 
function (|20)) represents a fully polarised spin-valley fer- 
romagnet. 

Similarly to Laughlin's wave function, the theory of 
composite fermions (CF)i^ can be extended to include 
the SU(4) internal degree of freedom.— Still, the main 
physical consequences of this additional symmetry can 
be captured within the simpler framework of Halperin 
wave functions4i Such states have been introduced in 
1983, soon after Laughlin's original work, to take into 
account the electronic spin and to describe non-fully po- 
larized FQHE states. This set of wave functions is readily 
generalised to the four-component case in graphene^^ 

(22) 

in terms of the product 

j=l kj<lj 

(23) 

of Laughlin wave functions for the four spin- valley com- 
ponents and the term 

4 Ni Nj 

cr=nnn(-i:^--gO"^' (^4) 

i<j ki kj 

which describes inter-component correlations. Here, z) 
is the complex coordinate of a particle in the compo- 
nent j, and Nj is the total number of j-type particles. 
As in the case of Laughlin's wave function, the power- 
counting argument relates the exponents rrij and to 
the component filling factors Vj = Nj/Nb- Indeed, for an 
arbitrarily chosen component j, the maximal exponent is 

Nb-1^ mj{Nj - 1) + n,jN,^j. (25) 

One notices that the inter-component correlations induce 
additional zeros in the wave function; this is energetically 
favourable because the SU(4)-symmetric Coulomb inter- 
action is as strong between particles of the same compo- 
nent as between those belonging to different ones. Fur- 
thermore, one notices that Eq. (P5|) has the character 
of a matrix equation, and it turns out to be useful to 
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introduce the exponent matrix = n^- the diagonal 
elements of which are simply the intra-component expo- 
nents Ujj = nij and the off-diagonal ones those corre- 
sponding to inter-component correlations. In terms of 
this exponent matrix, the relation between the compo- 
nent filling factors and the exponents reads 













1 






1 









In the zero-energy LL, the total filling factor is related 
to the component filling factors by 

4 

i/=-2 + ^z.„ (27) 

as a consequence of its half-filling for v =- whereas in 
all other LLs n the filling factor reads 

4 

i^ = ±[4(n-l) + 2]+^i/j. (28) 
i 

Notice that a state at a filling factor v is related to an- 
other one at — by particle- hole symmetry. 

In addition to the determination of the component fill- 
ing factors, the exponent matrix M. is useful also in 
two other respects. First, it allows one to distiguish 
between potential physical states and those that can- 
not describe a homogeneous liquid state that displays 
the FQHE. Indeed, the matrix must be positive defi- 
nite, i.e. contain only positive (or zero) eigenvalues, un- 
less the corresponding state is unstable and undergoes 
a phase separation between the different components4^ 
Second, The matrix M. encodes prominent properties 
of the quasiparticle excitations, such as their fractional 
charge and their statistics.— Finally, the rank of the ma- 
trix M encodes the SU(4)-fcrromagnctic properties of 
the different statcSfi^, In order to illustrate this point, 
we first mention that Eq. (|26|) is only well-defined if 
M. is invertible (of rank 4). This means that all com- 
ponent filling factors are fixed, and thus also all polari- 
sations which are simple combinations of these factors; 
e.g. the spin polarisation (in the z-direction) is sim- 
ply given by = {N/2){v^ji + v^^k' - n.^r - ^i^K'), 
whereas the valley-pseudospin polarisation reads Pz = 
{N /2){v^^K — v^.K' + ^i,K — i^i,K')- As an example, one 
may invoke the state with rrij = 3 for all j and = 2 
for all i ^ All component filling factors are fixed to 
be Vj = 1/9, as may be seen from Eq. p6)) . and this 
state would be an SU(4)-singlet candidate for a (yet un- 
observed) FQHE at = -2 + 4/9. 

In the opposite limit, where M. is of rank 1, the com- 
ponent filling factors are fully undetermined - the only 
combination that is fixed is the total sum. This is pre- 
cisely the case of Laughlin's wave function, which may 



be described as a four-component Halperin wave func- 
tion (|22p with all = mj = m being the same odd 
integer. As we have already mentioned, this corresponds 
to a fully polarised SU(4)-ferromagnetic state, and the 
sum of component filling factors is just Vj = 1/m. 

There are intermediate states for which, e.g., only one 
of the polarisations is fixed. As an example, we may con- 
sider the state with ruj = 3, ni2 = ni4 = 7123 = 7134 — 3, 
and ?ii3 = 7124 = 2, which is described by an expo- 
nent matrix of rank 2. In addition to the total filling 
factor, which is fixed sX v = —2 -f 2/5, the combina- 
tions 1^1 + 1^3 = 1/5 and z/2 + 1^4 = 1/5 are fixed. If 
we identify, for illustration reasons, the components as 
{1,2,3,4} - {{^,K),{^,K'),{UK),{i,K')Y this state 
would correspond to a valley-pseudospin singlet, with 

VK = i^(t,if) + i^(;,if) = 1/5, >^K' = ^(■\,K-) + >^(i,K') = 1/5, 

such that Pz = {N/2){vk — vk') = 0, whereas the spin 
is polarised and free to be oriented, e.g. by an external 
Zceman effect. 

Interestingly, the 2/5 and 4/9 states discussed above 
are in competition with completely polarised CF states, 
that may occur at fiUing factors = — 2 -I- p/{2p + 1), 
in terms of the integers p. Notice that the Halperin-type 
states at 2/5 and 4/9, which we discuss here, may alter- 
natively be viewed as unpolarised CF states.— Numeri- 
cal calculations have shown that the polarised states are 
generally higher in energy in the zero-energy graphene 
LLfi2, but they may become competitive when exter- 
nal symmetry-breaking terms are taken into account, 
such that one may expect similar spin transistions at 
fixed filling factors as in 2D electron systems in GaAs 
hetcrostructuresJ^ Quite generally, it is important to 
stress that the polarisation may change drastically when 
varying the filling factor - even if the LL degeneracy may 
be lifted in a precise hierarchy at the integer filling fac- 
tors 1^ = and h' = ±1, this hierarchy is easily destroyed 
when shifting the filling factor away from these values, 
as a consequence of the dominant Coulomb interaction. 
From a theoretical point one therefore expects that a fully 
spin- valley polarised state at i^=— 2 + 1/3 is depolarised 
when increasing the filling factor - this depolarisation 
is very efficient because of the low-energy skyrmion ex- 
citations of the SU(4) ferromagnetism associated with 
the Laughlin state. Indeed, upon increase of v one ob- 
tains a state at i' ~ —2 + 2/5 that is polarised in only 
one of the channels, e.g the spin for a Zeeman effect in 
the absence of vallcy-symmctry breaking terms. Upon 
further increase, the spin- valley polarisation disappears 
completely at ly = —2 + 4/9, where an SU(4) singlet is 
the state of lowest energy, in the absence of extrinsic 
symmetry-breaking effects. However, when approaching 
the filling factor v = —1, the SU(4) spin- valley polarisa- 
tion is again expected to be fully restored. 

We emphasise that, apart from the above-mentioned 
states at ^ = 2/5 and 4/9, most of the Halperin states 
(22) are not eigenstates of the SU(4) symmetry group and 
thus not of the bare SU(4)-symmetric Coulomb interac- 
tion. However, as we discuss in the following subsection 
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(a) K=l 




(c) k:=3 



-• « m- 




FIG. 2: Sketch of different states 4'v=—2+k,+i/3 of the 1/3 family, 
for (a) K = 1, (b) ft = 2, and (c) k = 3. 

for some of such states, they can be stabihsed with the 
help of relatively weak symmetry-breaking terms. We fi- 
nally notice that this picture and the expected polarisa- 
tion of the discussed FQHE states also holds true, within 
numerical calculations, fo the LL n = 1.— 




FIG. 3: Energy spectrum obtained in Ref. [sj by exact diagonal- 
isation of = 17 electrons interacting via the Coulomb interac- 
tion on a sphere threaded by Ng = 6 flux quanta. The spectrum 
is shown as a function of the Zeeman effect, in units of e^/dg. 
The colours indicate levels with the different spin polarisations 
52 = 11/2 (green, corresponding to the state V1/3) 9/2 (blue) 
adn 7/2 (pink). Rod dot correspond to energy levels with other 
polarisation. 



1. The 1/3 family of FQHE states 

We have discussed, until now, the use of Halperin wave 
functions for potential FQHE states in the filling-factor 
range —2 < v < —1. Notice, however, the the only 
n = FQHE states that have been clearly identified 
a.r ly ~ ±1/3, ±2/3, and ±4/3, whereas the states at 
v ~ ±5/3 are absent or extremely weak.— ^— The un- 
derstanding of these states requires the inclusion of fully 
occupied spin-valley subbranches in addition to partially 
filled ones. In principle, these states may also be de- 
scribed in terms of generalised Halperin wave functions 
with broken SU(4) symmetry, and we terminate this re- 
view with a brief discussion of them. The states, which 
are depicted in Fig. [2l may be constructed by the wave 
functions 

K Nb Nb/3 

^.=-2+.+i/3 = n n (4f n k-^o\ 

j kj<lj k<l 

(29) 

(i) 

where are the complex coordinates of j-type particles 
that reside in the k fully occupied spin-valley branches 
J = !,...,/« (since Vj = 1, we have Nj = Nb), whereas 
Wk is that of a particle in the other spin-valley branches, 
occupied by Nb/S particles, where we do not specify ex- 
plicitly the component. This state is described by an 
exponent matrix with mj ~ 1 for j — uij = 3 

for j = K + 1, 4, riij = if one of the indices is 1, k 
and Hij = 3 otherwise. In the framework of quantum 
Hall ferromagnetism, the state may be viewed as k "in- 
ert" levels,— whereas the electrons in the 4 — k partially 
occupied subbranches form a Laughlin-type state with 
incorporated SU(4 — k) ferromagnetic low-energy exci- 
tations in terms of (3 — K)-fold degenerated spin-valley 



waves. Other members of the 1 /3 family may be obtained 
from the states ((29| with the help of a particle-hole trans- 
formation. 

As discussed in the previous section, the states 
cannot describe the ground state of the Coulomb in- 
teraction because they are no eigenstates of the SU(4) 
symmetry. These states may be stabilised artificially by 
a particular choice of the interaction potential between 
the different particles that explicitly breaks the SU(4) 
symmetry.— However, the more physical approach which 
we adopt here shows that the states ([29|) may be rele- 
vant even for an SU(4)-symmetric interaction potential 
as soon as the extrinsic symmetry-breaking terms are in- 
cluded. 

For K = 2 [see Fig. HKb)], i.e. for a filling factor 
ly = 1/3, the energy spectrum obtained from exact di- 
agonalisation is shown in Fig. [3] as a function of an ex- 
trinsic symmetry-breaking field that has been chosen to 
be the Zeeman effect. The spectrum has been obtained 
with the help of the DiagHaniS code with implemented 
SU(4) symmetry for up to = 17 electrons on a sphere 
threaded by Nb = 6 flux quanta. As expected from the 
above discussion, the ground state in the absence of a 
Zeeman effect is not the state ■(/'i/=i/3 described by Eq. 
([29| because it does not have the correct spin polarisa- 
tion. However, the state is stabilised already for very 
small symmetry breaking, i.e. for a Zeeman effect above 
— O.Ole^/eZs, as may be seen in Fig. [31 where the 
ground state has the correct (maximal) spin polarisation 
Sz = 11/2 (green dots). Although the state ?/'i/=i/3 is the 
ground state above this critical value of the Zeeman ef- 
fect, its low-energy excitations are not the usual collective 
excitations of the Laughlin state, but coherent spin-flip 
excitations in the sector Sz =9/2 that are represented 
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by the blue line in Fig. |3l These spin-flip excatiations are 
the relevant modes below a second critical Zeeman field 
~ O.OSe^/eZs, whereas above A| the lowest-energy 
excitations are the usual charge excitations in the same 
polarisation sector. 

The case k = 1, which corresponds to a single fully 
occupied spin- valley branch [see Fig. Ula)], may be 
checked within a simplified two-component scheme that 
neglects the two-fold degenerate spin-valley modes asso- 
ciated with the partially filled three subbranches j = 2, 
3 and 4. The corresponding (simplified) wave function 
reads 

Nb Nb/3 

V'^=-°iTi/3 = n (^'^ ~ n -wif^ (30) 

k<l k<l 

where Zk is the position of a particle in the fully occu- 
pied subbranch and Wk one in the 1/3- filled second com- 
ponent. This wave function has been tested in exact- 
diagonalisation calculations with the help of the Di- 
agHam code^ with an implemented SU(2)-symmetric 
Coulomb potential, for N = 22 electrons on a sphere 
with Nb = 15 fiux quantai^ The obtained energy spec- 
trum shows the same features as that depicted in Fig. 
[3] obtained within a four-component calculation, with a 
similar critical field ~ O.Ole^/e/s above which the 
state is stabilised albeit with a slightly larger field 
A^ ~ O.OSe^/e/s, below which the collective excitations 
are dominated by spin-fiip excitations.— 

These two results obtained numerically, for n = 2 and 
for K = 1 in a simplified version, hint at a certain univer- 
sality in the mechanism of stabilising states of the form 
(|29p by weak extrinsic symmetry-breaking fields. The 
physical picture that emerges from it may be summarised 
as follows: the systems has an interaction-driven ten- 
dency to form such states (here the states ip^^_2+K,+i/3 
of the 1/3 family), but a small extrinsic SU(4) spin- valley 
symmetry breaking is nevertheless necessary to stabilise 
them. This needs to be contrasted to the SU(4) quan- 
tum Hall ferromagnetism discussed in Sec. IIII Al where 
the state remains stable even in the complete absence of 
extrinsic symmetry-breaking effects. 

We finally notice that, even in the intermediate regime 
A^ < Az < A|-, the lowest-energy excitations in the 
limit qIb ^ ^ arc not the collective spin- flip excitations, 
but as for the usual Laughlin 1/3 state quasi-particle 
excitations that may eventually be dressed by SU(4 — k) 
spin- valley textures in the partially occupied subbranches 
and that are responsible for the activation gap measured 
in the experiments J^iii 

IV. CONCLUSIONS 

In conclusion, we have reviewed theoretically the role 
of electronic interactions in graphene Landau levels. 



These interactions are responsible for two prominent ef- 
fects: (a) the formation of SU(4) spin-valley quantum 
Hall ferromagnets that are likely to be responsible for 
the observed IQHE at z/ = 0, ±1, ±3, ±4, and ±5 that 
do not belong to the series ([T]) of the usual (relativistic) 
graphene IQHE; and (b) the recently observed FQHE. 
Although both effects arc known also in the context of 
non-relativistic quantum Hall systems, such as in GaAs 
heterostructures, they are different in graphene as a con- 
sequence of the approximate SU(4) spin- valley symmetry 
of the Coulomb interaction potential. Even if the SU(4) 
symmetry of graphene LLs is broken by extrinsic effects, 
such as the Zeeman effect or a valley-pseudospin Zeeman- 
type effect due to static lattice distortions in graphene, 
the latter effects are associated with energy scales that 
are much smaller than the leading Coulomb interaction 
scale, for physically accessible magnetic fields. These ex- 
trinsic effects are mainly cooperative with the tendency 
of forming maximally spin-valley polarised states. In 
the context of quantum Hall ferromagnetism, they orient 
the preformed spin-valley magnetisation into particular 
channels, whereas they are necessary to stabilise the trial 
states 'ipi^^-2+K+i/3 that may account for the experimen- 
tally observed members of the 1 /3 family. 

Also other FQHE states than those of the above- 
mentioned 1/3 family may be described in the framework 
of the SU(4) theory of the FQHE and are expected to dis- 
play very special spin- valley polarisations. It remains an 
experimental challenge to have access to these states and 
their physical properties, but from an experimental point 
of view we seem to be only at the beginning of the dis- 
covery of the possibly very rich physical properties of the 
graphene FQHE. The expected findings of novel FQHE 
states in graphene may provide other surprises that will 
certainly also challenge the SU(4) theory of the graphene 
FQHE. 
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